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The volume and surface symmetry parts of the nuclear symmetry en- 
ergy and other coefficients of the liquid droplet model are determined from 
the measured atomic masses by the maximum likelihood estimator. The 
volume symmetry energy coefficient extracted from finite nuclei provides 
a constraint on the nuclear symmetry energy. This approach also yields 
the neutron skin of a finite nucleus through its relationship with the vol- 
ume and surface symmetry terms and the Coulomb energy coefficient. The 
description of nuclear matter from the isoscalar and isovector components 
of the density dependent M3Y effective interaction provide a value of the 
symmetry energy that is consistent with the empirical value of the symme- 
try energy extracted from measured atomic masses and with other modern 
theoretical descriptions of nuclear matter. 
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1. Introduction 

The investigation of constraints of nuclear symmetry energy has recently 
received new impetus with the plans to construct a new accelerator facility 
(FAIR) at GSI Dramstadt. The nuclear symmetry energy (NSE) plays a 
central role in determining to a large extent the equation of state (EoS) 
of isospin asymmetric nuclear matter, proton fraction in neutron stars and 
neutron skin in heavy nuclei and it enters as an input to the heavy ion 
reactions |Tj , [2] . Various many body calculations using a realistic nucleon- 
nucleon (NN) interaction as input (Brueckner-Hartree-Fock [3] or Dirac- 
Brueckner-Hartree-Fock [3] and the variational method [5] etc.) lead to 
rather different results for the symmetry energy. In view of the large dif- 
ferences present between various calculations of the symmetry energy even 
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at subsaturation densities, the question arises naturally whether one can 
obtain empirical constraints from finite nuclei. As the degree of isospin 
diffusion in heavy-ion collisions at intermediate energies is affected by the 
stiffness of the nuclear symmetry energy, these reactions provide constraints 
on the behaviour of the nuclear symmetry energy at subsaturation densi- 
ties [6J . Traditionally, the symmetry energy of nuclei at saturation density 
is extracted by fitting ground state masses with various versions of liquid 
drop mass formula (LDM). To this end one needs to decompose the sym- 
metry term of LDM into the bulk (volume) and surface terms [7] along the 
lines of the liquid droplet model and identify the volume symmetry energy 
coefficient as the symmetry energy derived from infinite nuclear matter at 
saturation density. The coefficients of liquid droplet model are extracted by 
employing the maximum likelihood estimator method. 

In the present work, the nuclear symmetry energy is calculated theo- 
retically using the isoscalar and isovector components of M3Y-Reid-Elliott 
effective interaction supplemented by a zero range pseudo-potential along 
with the density dependence (DDM3Y) and its value at saturation density 
is compared with the volume symmetry energy coefficient extracted from a 
fit to the atomic mass excesses from the latest mass table [8]. The M3Y 
interaction was derived [9J by fitting its matrix elements in an oscillator 
basis to those elements of the G-matrix obtained with the Reid-Elliott soft- 
core NN interaction. The ranges of the M3Y forces were chosen to ensure a 
long-range tail of the one-pion exchange potential as well as a short range 
repulsive part simulating the exchange of heavier mesons [ID]. The zero- 
range pseudo-potential represented the single-nucleon exchange term while 
the density dependence accounted for the higher order exchange effects and 
the Pauli blocking effects. The real part of the proton- nucleus interaction 
potential obtained by folding in the density distribution function of the in- 
teracting nucleus with the DDM3Y effective interaction is found to provide 
good descriptions of elastic and inelastic scatterings of high energy protons 

and proton radioactivity [12] . In the present work, the results obtained 
for the symmetry energy from DDM3Y interaction are consistent with the 
results obtained by the fitting of masses to the liquid droplet model. 

2. The nuclear symmetry energy 

The nuclear EoS can be expanded in terms of isospin asymmetry X as 

e(p, X) = e(p, 0) + S(p)X 2 + S l {p)X A + S 2 (p)X (i + (1) 

which provides the standard definition of the nuclear symmetry energy S(p) 
where e(p, X) is the energy per nucleon of nuclear matter with isospin asym- 
metry X = P p " + P p p p , p = Pn+Pp, where p n , p p and p are the neutron, proton 
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and nucleonic densities respectively. The DDM3Y effective NN interaction 
is given by 

v 00 (s,p,e) =C 3Y (s,e)g(p), v i(s, p, e) = t^ Y (s, e)g(p) (2) 

where s is distance between two interacting nucleons, g(p) = C(l — /3p 2//3 ) 
is the density dependence and the isoscalar tQQ 3Y and the isovector tg( 3y 
components of M3Y interaction potentials [ID], [T2] supplemented by zero 
range potentials are given by 



t^ Y (s, e) = 7999^^1 - 2134^^ - 276(1 - ae)5(s) (3) 
4s 2.5s 



and 



iff* (s, e) = -4886 6Xp( ~ 4g) + 1176 ex P( 2 ' 5g ) + 228(1 _ ae)5(s) (4) 

4s 2.5s 

respectively, where the energy dependence parameter a=0.005/MeV. Based 
on the Hartree or mean field assumption and using the DDM3Y interaction, 
the EoS can be derived as 

e( P ,X) = l^MX) + (^)(1 " Pp 2/ % (5) 

where k F = 1.57r 2 p, J v = J v00 + X 2 J v01 , J v00 (e) = J J J t^ 3Y (s, e)d 3 s and 
JvOl( e ) = 1 1 1 t()[ 3Y (s,e)d 3 s represent the volume integrals of the isoscalar 
and isovector parts of the M3Y interaction and 

F(X) _ p+xroni-^ (6) 

where m is the nucleonic mass equal to 938.91897 MeV/c 2 . In nuclear 
matter, e(p,X) = e(p,0) + e'(p,X) so that e'(p,X) = e(p,X) - e(p,0) sa 
S(p)X 2 for other terms being small. Therefore, an alternative physical 
definition of the nuclear symmetry energy [13] is the energy required per 
nucleon to change the symmetric nuclear matter (SNM) to pure neutron 
matter (PNM) which is given by 

S(p) = e PNM - esNM = (2 2/3 - 1)^(-^) 2 / 3 + %p(l - (3p 2 l 3 )J v01 (7) 

5 po 2 



where cpnm = e(p, 1), £SNM = e (p>0) are the energy per particle for PNM 
and SNM respectively, po is the saturation nucleonic density, E F 



h 2 k% 

2m 
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is the Fermi energy for the SNM in the ground state with /cp as the cor- 
responding Fermi momentum. The constants of density dependence C and 
(3 of the effective interaction are obtained by reproducing the saturation 
energy per nucleon and the saturation density of SNM [TJj. In ref. |14j the 
symmetric nuclear matter properties such as its EoS and incompressibility 
were calculated whereas in this work we have extended that for asymmetric 
nuclear matter in order to calculate the nuclear symmetry energy. 

The first term of the right hand side of Eq.(7) is the kinetic energy 
contribution whereas the second term is the potential energy contribution 
and accounts for the nuclear interaction. If one uses definition of Eq.(l) 
for the nuclear symmetry energy, then the second term remains unaltered 
while the the first term reduces by about five percent to Ep(-j^) 2 / 3 ]. 

The standard definition of Eq.(l) also provides Si(p) = Jr[f #f(^) 2/3 ] 
and S-2,(p) = §§[§ Ep(-j^) 2 / 3 ] etc. Interestingly, the definition of nuclear 
symmetry energy e sym given in ref. [15] yields S(p) = 0.5e syrn = 28 MeV and 
Si(p) = —S(p)X = -18.76 MeV which are independent of nuceonic density. 



3. The liquid droplet model of nuclei and symmetry energy 

The volume and surface terms in the standard semi empirical mass for- 
mula pertain to the isospin symmetric systems. The volume coefficient 
provides the binding energy per nucleon whereas the surface coefficient, up 
to a certain extent, provides the surface energy. The symmetry term in 
the standard binding energy formula has a volume character only. But 
when the surface energy is affected by the isospin asymmetry, the ther- 
modynamic consistency requires that some of the asymmetry moves to the 
surface. Minimization of the net nuclear energy with respect to the par- 
titioning of asymmetry produces an expression [7j for the binding energy 
B(A,Z) of a nucleus with mass number A and the atomic number Z given 

by 



B(A, Z) — a v A - - a c ^> - ; - Z± > + 6, (8) 

Ss 

where 5 = a p A~ l l 2 for even N-even Z, — a p A~ l l 2 for odd N-odd Z, for 
odd A, and neutron number N=A-Z. The above expression is similar to the 
droplet model where skin size is a basic parameter and one of the starting 
points. S v and S s are now the volume and surface symmetry parameters, re- 
spectively, whereas a v , a s , a c and a p are the usual volume, surface, coulomb 
and pairing energy coefficients. Allowing the mass number A going to in- 
finity, it may be seen that the volume symmetry energy coefficient S v is 
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equal to the NSE obtained from the (infinite) nuclear matter calculation. 
Therefore extracting S v from measured atomic mass excesses provide ex- 
perimental value for the NSE at normal nuclear density. Theoretical atomic 
mass excesses AMaz can be obtained from the theoretical binding energy 
B(A,Z) by correcting for the electronic binding energy as 



AMa,z = ZAm H + {A- Z)Am n - a el Z 2 - 39 - b d Z™ 5 - B(A, Z) (9) 

where Amu = m p + m e — u = 7.28897050 MeV + a e i + b e i and Am n = 
m n — u = 8.07131710 MeV, m p , m n , m e are the masses of proton, neutron 
and electron and u is the atomic mass unit, all expressed in MeV and the 
electronic binding energy constants [16] a e i = 1.44381 x 10~ 5 MeV and 
b e i = 1.55468 x 10 -12 MeV. This approach [7] also yields a relationship 
among neutron skin, a c , S v and S s given by 

Rn-R P _ A N-Z-^ZA^ 

R 6NZ 1 + ^AV3 1 ' 

The difference between equivalent sharp radii for neutrons R n and protons 
R p is primarily linear in the asymmetry and the symmetry coefficient ratio 
S v /S s measures the neutron skin of a nucleus. 



4. Results and discussion 

We extract the nuclear symmetry energy S v from measured atomic mass 
excesses and associated errors using the maximum likelihood method de- 
scribed in detail in ref.[T7]. This leads to the generalised equations [17]. 



£ [ AM lx -(AM lh+ ^ ]d{ AM k)= ^ m 

t[ ai x 2 + a th ^ d Pu 

" [AMl x - {AMl h + y th *)f - (*l x 2 + a th 2 *) _ Q 

" [AMl x -{AMl h + ii th *)\ ^ (i3) 

i= i (vix 2 + <?th 2 *) 

where p u are the unknown m parameters of the model. Here AM\ X is the 
measured mass excess for a particular nucleus for proton number Z and 
neutron number N, and AAfi is the corresponding calculated quantity 
and a* is the associated error in each of n such measurements, ath is 
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the intrinsic model error which accounts for known and unknown missing 
terms in the theoretical model used for fitting the mass excesses. Here we 
assume that the true mass excess u\ r of the nucleus i can be written as 
u\ r = AM£ h + e\ h , where e\ h is the theoretical error term and is distributed 
normally as e\ h € N(fj, t h, a t h) with a mean fi t h and a standard deviation 
o t h around this mean. The notations a t h 2 * and fi t h* mean that by solving 
Eqs.(12,13) we obtain the estimates ath 2 * an d fith* of the true errand [ith- 
Use of root-mean-square deviation (a rms ) defined as 



1 n 

-V(am;-a^) 2 ] 1/2 (i4) 



n 

II ^ 

t=l 



as the error of the theoretical mass model and obtained by minimising its 
value by adjusting the model parameters is reasonable when all the errors 
a\ x associated with the measurements are small compared to the model error 
a rms . However, for large experimental errors a\ x this definition is unsatis- 
factory, since both the theoretical and the experimental errors contribute 
to the rms deviation. We must therefore use an approach that decouples 
the theoretical and the experimental errors. In the present case, since it is 
a minimally modified formula without the shell corrections or the Wigner 
term, the theoretical model error must account for the various other known 
and unknown terms in the model. The model error (ath) obtained in this 
way (ML estimation) contains no contributions from the experimental un- 
certainties a l ex . Thus, we have two additional equations here compared to 
usual least square equations (minimising a 2 or alternatively x 2 ) that arise 
when model parameters are estimated by adjustments to experimental data 
under the assumption of a perfect theory with ath = and fj,th = 0- The 
above equations are equivalent to minimizing S with respect to p v , where 



* = E 



i=i 



°lx + a th 



2* 



and solving 



where 



1 n 

oth 2 * = — ~ a E PC - A< - n t h*f - al x 2 } (16) 

Ei=i^ i=i 

»th* = — - jr E u£" [(AMi x - AO] (17) 

E"=i^i i=i 



(18) 



(ai x z + a th 2 *) k 
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K = 2 
k„ = 1 



(19) 
(20) 



The unknowns nth* and ath 2 * are then determined from Eqs.(16,17) by an 
iterative procedure whose convergence is found to be quite good. This way 
the experimental error is subtracted from the difference between the exper- 
imental and the calculated mass excesses. As the model considered here 
does not contain any term like clqA , that is strictly a constant parameter, 
the most complete characterisation of the theoretical error requires both its 
mean [x t h and its standard deviation ath around this mean. Hence we need 
to solve the full m + 2 set of equations. If ^ t h* is found to be significantly 
different from zero the theory will need modification. 

In ref. [14] both the chi-square and the sum of deviation squares were 
minimized. The results of those two minimization did differ slightly. How- 
ever, the data are not expected to approach the theory when measurement 
errors tend to zero. In the present work the theoretical errors are assumed 
to accompany the experimental errors. In fitting experimental data, the 
theoretical errors are estimated simultaneously with the optimal parameter 
values. This procedure effectively produces a minimization which is inter- 
mediate between the minimization of the sum of deviation squares and that 
of the chi-square. This procedure gives the possibility for a realistic estima- 
tion of the uncertainties in the fitted parameters. As one may compare with 
ref. [13] to find that the volume symmetry energy coefficient obtained in the 
present work is substantially different. This is due to the fact that Bethe- 
Weizsacker mass formula is minimally modified along the lines of the liquid 
droplet model by partitioning the symmetry term into volume and surface 
terms. Other energy coefficients do not differ much beyond the limits of the 
earlier work [14j. 

The coefficients of the liquid droplet model [Eq.(8)] are evaluated by 
fitting the recent measured and extrapolated atomic mass excesses from 
Audi-Wapstra-Thibault atomic mass table [8] by minimizing S of Eq.(15) 
and are shown in Table- 1. The ath and nth simultaneously solved from 
Eqs.(16,17) are also tabulated. The values of a t h and nth f° r 2228 experi- 
mentally measured atomic mass excesses are 2.880 and 0.029 respectively. 
Exclusion of the measured atomic mass excesses of lighter nuclei having mass 
number A < 16 results in ath = 2.782 and fith = 0.037. When the additional 
951 extrapolated data are included for the same analysis, that is for total 
3179 measured+extrapolated data, the values obtained are ath = 2.960 and 
fJ-th = 0.040. These values are acceptable as the values of fith do not differ 
significantly from zero. Fig.l and Fig. 2 show the plots of fitting errors of 
atomic mass excesses versus mass numbers for the droplet model [Eqs.(8,9)] 
mass formula. 
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Table 1. Coefficients of the liquid droplet model mass formula extracted from 
atomic mass excesses. 



a v 
MeV 


a s 
MeV 


a c 
MeV 


S v 
MeV 


S s 
MeV 


a p 
MeV 


[a] 15.563 
±0.00043 


17.652 
±0.00107 


0.695 
±0.00004 


29.687 
±0.0097 


17.680 
±0.0216 


11.813 
±0.0044 


[b] 15.500 
±0.00012 


17.480 
±0.00025 


0.689 
±0.00002 


30.048 
±0.0042 


16.674 
±0.0071 


10.246 
±0.00048 


[c] 15.465 
±0.00012 


17.394 
±0.00025 


0.686 
±0.00002 


30.130 
±0.0043 


16.317 
±0.0067 


10.273 
±0.00049 



[a]. Using experimentally measured 2168 atomic mass excesses for A>16. 
p th = 0.037 and a th = 2.782 
[b]. Using all the experimentally measured 2228 atomic mass excesses. 
Hth = 0.029 and a th = 2.880 
[c]. Using measured 2228 + extrapolated 951 atomic mass excesses. 
Pth = 0.040 and a th = 2.960 

A mean field calculations with DDM3Y effective interaction is performed 
using the usual values of energy dependence parameter a = 0.005MeU _1 
[TP] , the saturation density po = 0.1533/m -3 [T2] and the saturation energy 
per nucleon eo = —15.26 ± 0.52 MeV. This saturation energy per nucleon is 
the volume energy coefficient and this value used in the present calculations 
more or less covers the entire range of values of a v obtained here (Table- 1) 
or fitting the original Bethe-Weizsacker mass formula [14J or other droplet 
models [13 13 [201 [2H [22] . The value obtained for the NSE at the saturation 
density, S(po), is found to be 31.18 ± 0.29 MeV using definition of Eq.(l) 
and 31.86 ± 0.29 MeV using definition of Eq.(7). 

The value of S v = 30.048 ±0.004 MeV extracted from experimental mass 
excesses is reasonably close to the theoretical estimate of the value of NSE 
at the saturation density S(po) described above. For A > 16 little change 
is observed to the value of S v which becomes 29.687 ± 0.010 whereas the S s 
changes to a larger extent to 17.680 ± 0.022. This behaviour suggests the 
fact that the surface energy depends upon symmetry. The value obtained 
for S v in ref. [7] is between 29.10 MeV to 32.67 MeV and that obtained by 
the liquid droplet model calculation of ref. [19] is 27.3 MeV whereas in ref. 
|23j it is 28.0 MeV. It should be mentioned that the value of the volume 
symmetry energy co-efficient S v in some advanced mass description |24] is 
close to the present value which with their —K vo i.b vo i = S v equals 29.3 MeV. 
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Fig. 1. The plot of differences between 2168 (A>16) measured and theoretical 
atomic mass excesses calculated by the liquid droplet model [Eq.(8)] mass formula 
versus mass number A. 



In a very recent similar work [25] the volume symmetry energy co-efficient 
comes out to be about 29 MeV when with other terms, pairing energy, 
Wigner term and shell corrections are also included. The ratio S v / S s which 
is a measure of the neutron skin thickness is found to be about 1.8 in the 
present calculations. The value of this ratio obtained in ref. [7J is about 
2.0 to 2.8 whereas in liquid droplet model calculations of ref. [19] the value 
obtained is 1.68 and in ref. [23] it is calculated to be 1.3. 

The value of NSE at nuclear saturation density ~ 30 MeV, therefore, 
seems well established empirically. However, theoretically different sets of 
parametrizations of the relativistic mean-field (RMF) models, which fit the 
obseravables for isospin symmetric nuclei well, lead to a relatively wide range 
of predictions 24-40 MeV for S(po). In Table-2 results for the S(po) using 
DDM3Y interaction are compared with the results from the variational cal- 
culations using the Argonne and Urbana NN potentials, in combination 
with Urbana models for the three- nucleon interaction [TAT]. The last col- 
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Fig. 2. The plot of differences between measured 2228 and theoretical atomic mass 
excesses calculated by the liquid droplet model [Eq.(8)] mass formula versus mass 
number A. 

Table 2. Result for the S(po) of the present mean field calculation is compared 
with the results from the variational calculations of Rcfs. [5, 26 using the Argonne 
and Urbana NN potentials, in combination with Urbana models for the TNI. The 
last column includes a relativistic boost correction Sv and the adjusted UIX* TNI. 



Present Calc. 


Avl4 


Avl4+UVIII 


Uvl4 


MeV 


MeV 


MeV 


MeV 


[a] 31.18 ±0.29 


24.90 


27.49 


26.39 


[b] 31.86 ± 0.29 




Uvl4+UVIII 


Avl8 


Avl8+UIX 


Avl8+<fa+UIX* 


MeV 


MeV 


MeV 


MeV 


28.76 


26.92 


29.23 


30.1 



[a] . Using definition of Eq.(l) for the nuclear symmetry energy. 

[b] . Using definition of Eq.(7) for the nuclear symmetry energy. 
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umn includes a relativistic boost correction and the adjusted UIX* TNI. 
The present result of the mean field calculation is close to the result of 
Avl8+<k>+ULX* variational calculation [5]. 

5. Summary and conclusion 

In summary, we show that theoretical description of nuclear matter 
based on mean field calculation using density dependent M3Y effective NN 
interaction gives a value of the symmetry energy that is consistent with the 
empirical value extracted by fitting the droplet model to the measured and 
extrapolated atomic mass excesses using the maximum likelihood estima- 
tor method. The value of volume symmetry energy coefficiet changes little 
compared to the surface symmetry energy coefficient when measured atomic 
mass excesses of nuclei lighter than 16 are excluded. This observation high- 
lights the fact that the symmetry energy depends on the surface tension 
and vice versa. The volume and surface symmetry energy coefficients are 
related to the neutron skin of finite nucleus. Such mean field calculations of 
nuclear symmetry energy thus satisfy the constraints from finite nuclei and 
also agree with recent theoretical descriptions of nuclear matter. 
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